Electrodynamics of a split-ring Josephson resonator in a microwave line 
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We consider the coupling of an electromagnetic wave to a split-ring Josephson oscillator or radio- 
frequency SQUID in the hysteretic regime. This device is similar to an atomic system in that it 
has a number of steady states. We show that one can switch between these with a suitable short 
external microwave pulse. The steady states can be characterized by their resonant lines which 
are of the Fano type. Using a static magnetic field, we can shift these spectral lines and lift their 
degeneracy. 

PACS numbers: Josephson devices, 85.25.Cp, Metamaterials 81.05.Xj, Microwave radiation receivers and 
detectors, 07.57.Kp 



I. INTRODUCTION 

In the past decade artificial materials (metamaterials) 
have been developed by inserting metal split-rings or rods 
as metaatoms into dielectric materials. This way nega- 
tive index materials have been fabricated Q, @- These 
metaatoms spaced periodically enable the medium to 
have an engineered electromagnetic response. In prin- 
ciple such an artificial atom has many advantages over 
a real atom. In real materials, the dipole moment of 
the atom is very small and the interaction with electro- 
magnetic field is small. Therefore to realize significant 
interactions, large densities are necessary. On the con- 
trary an artificial atom can be constructed with a high 
dielectric or a high magnetic dipole moment so that it 
will respond strongly to electromagnetic waves. 

Many candidates for artificial atoms have been pro- 
posed, most of them with intrinsic nonlinearities. Among 
them are resonant conductive elements with inserted 
(nonlinear) diodes Q, split-ring resonators loaded with 
varactor diodes 0], Kerr materials @ or laser amplifiers 
Q. Unfortunately many of these systems exhibit large 
losses. To remedy this problem it has been been sug- 
gested to work with superconducting split-ring resonator 
arrays. An example is the study [7| using from high- 
temperature superconducting films. See also the review 
8] for more details. 

Among these superconducting meta-materials are the 
split-ring Josephson resonators subsequently referred to 
as SRR-JJ where a Josephson junction, a superconduct- 



ing weak link, is embedded in the slit. In the Joseph- 
son community, this device is called an RF SQUID for 
Radio Fre quen cy Superconducting Quantum Interference 

Device HG3. 

The device we consider is shown in the left panel of 
Fig. [T] It is a split ring resonator in which is embedded 
a Josephson junction. Practically it can be made using a 
ring like strip of superconducting material where a small 
region was oxidized to make the junction. The right panel 
of Fig. [1] shows the electric representation of the device, 
an inductance L for the strip and the Resistive Shunted 
Junction (RSJ) model for the Josephson junction (JJ). 
The latter represents the junction as a resistor R, a ca- 
pacity C and the nonlinear element in parallel. This 
last element is the sine coupling I c sin(<I>/<I>o) where <E> is 
the magnetic flux and <E>o is the reduced flux quantum 
$ = ft/(2e). 

The use of these devices in metamaterials was advo- 
cated by Lazarides [ll| [l2j and by the authors [l3| . The 
electrodynamics of such metamaterials where split-rings 
are coupled inductively was investigated in [TJ [T^]- It 
leads to a 2D discrete sine-Gordon equation. For weak 
external radio- frequency forcing and small densities, this 
dilute gas of artificial atoms demonstrates a nonlinear 
magnetic response 

Such arrays of SRR-JJ arc difficult to analyze because 
the coupling between the artificial atoms is partly induc- 
tive, partly capacitive. Another point is that the field 
couples to the phase of the Josephson junctions so that 
the system should have two components. For these rea- 
sons, in this article we study the dilute system where 
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each SRR-JJ can be considered as isolated although cou- 
pled to an electromagnetic wave. We consider the so- 
called hysteretic regime where the system has controlled 
metastable states and show that one can switch from the 
ground state to one of these excited states by applying a 
suitable flux pulse. Then we can detect these states by 
examining the reflection coefficient of an electromagnetic 
wave incident on the device, this is a kind of spectroscopy. 
More precisely, the resonance observed are asymetric, of 
the Fano type. This is interesting because it is observed 
for many plasmonic systems [l5j. 

The article is organized as such. In section 2 we present 
the basic equations of the model and analyze it's steady 
states. In section 3 we show how to switch from one state 
to another by applying a suitable flux pulse. Section IV 
discusses the scattering of an electromagnetic wave by a 
SRR-JJ and this leads to a spectroscopy of the steady 
states. 



II. THE MODEL 

We consider here that the SRR-JJ is subject to irradi- 
ation by a microwave field. The equations describing the 
system light-ring are the generalized pendulum equation 
for the flux and the Maxwell equation for the electro- 
magnetic field. The Maxwell equations can be written 
as 

VxE = -B, t , V-D = 0, 
V x B = E 4) V-B = 0. 

where E and H are respectively the electric field and the 
magnetic field, and the subscritpts indicate the partial 
derivative. The fields B, H and E, D are related by 

B = i uoH + M, D = e E + P, 

where M and P are respectively the magnetization and 
polarisation of the medium. In the dispersionless limit we 
can write D = eoE~E . The magnetization M is located 
in the thin plane layer that is embedded into a dielectric 
surrounding with permittivity e at a point z — 0. The 
magnetization can then be represented as 

M(t,r) = m(t)n r l5(z), 

where I is the film thickness, n r is the density of SRR-JJ 
and m(t) is the magnetic moment of the ring ;7ii| . 

We assume that the plane electromagnetic wave prop- 
agates along the z-axis and is polarized so that H is par- 
allel to the y axis, the normal to the plane of the SRR-JJ 
and E is parallel to the x axis (see Fig. [TJ. Since all the 
magnetic moments of the SRR-JJ in the layer are paral- 
lel, the magnetization is parallel to H. In this case the 
Maxwell equations take the form 



They imply the wave equation 

H >zz - nos eH tU = e M^t (1) 
Ohm's law applied to the ring gives 
£ = Ul + U c , 

where £ is the electromotive force induced by the mag- 
netic field of the electromagnetic pulse incident on the 
ring. 

£ = -^«-S# (t (0,t), 

and where S is the surface enclosed by the ring. We 
neglect the resistance of this loop which we assume to 
be made of superconductive material. Following Joseph- 
son's second relation, the voltage across the junction is 

° 2e dt ° dt ' 

where tp is the phase of the junction. The inductance L 
due to the loop of the SRR results in the voltage 



U L =L 



dt 



where / is the current in the loop. Taking into account 
Ohm's law together with these definitions we obtain 

$ + § Q tp + LI = 0, 



or 



I = -L-\$ + $ <p). 



(2) 



The Resistive Shunted Junction model for the Joseph- 
son junction together with Kirchoff 's law for the current 
imply 



It 



c -t- *J 



C 



dUc 
dt 



+ I c sin ip + I R 



= C *°W +IoBm,p + ll8t 
Thus the variable tp evolves according to 



d 2 ip $ dtp T . 



-L 1 ($ + $ ^)- 



In this equation the flux $ has two components, the in- 
cident magnetic flux and the flux induced by the current 
in the loop of the SRR-JJ. That is 

$ = SH m (t) 



E 



-(HoH t + M it ), 



-EqEEj 



2 Y dt 
Since M — n r lSI(t) we have 

M = —n r lSL~ l (<b + $o(p). 
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The wave equation ([1} can be rewritten as 

where the one ring magnetization m can be defined as 

m(t) = n r lSI(t). 

Taking into account the definition of the current ([2} one 
can write 

rh(t) = -n r lSL-\<I> + <f> tp) = -rirlSL-^SH + $ ip). 
The phase variable ip is governed by the equation 

Ir . 1 SH(0,t) 




d 2 ip 1 d(p 
W + RC~dt 



• sin op 



C$ r LC $ 
The natural units of time, flux and space are given by 

c 



lj t = 1/vLC, $ , zo = 



where lot is the Thompson frequency and z is the inverse 
of the Thompson wave number. Hence, we can introduce 
the new dimensionless variables 

h = SH/$o, t = u/xt, Q = z/zq. 

In terms of these variables the final equations are 



h T7 



-n(h + v),tt$(C)) 



(3) 



(4) 



f,TT + Oiip. T + <p + f3 sin if = —h, 
where the parameters a, j3 and k are 

1 jL LI c n r lSeoUJT 

The system of equations (|3I4[) is our principal model 
and we will investigate it in detail throughout the article. 

Let us briefly recall the main properties of the SRR- JJ 
equation Q. It can be written as the 1st order system 



<f,r = "0, 

-atp — P sha(y) — ip + h, 



whose fixed points are (0, 0) and (ip* , 0) where 
- /3sm(tp*)-ip* + h = 0. 



(5) 
(6) 



(7) 



A plot of the above relation indicates that for f3 < 4.34, 
the equation has no solution so that there is only the 
(0, 0) fixed point. For large (3, the fixed points can be 
approximated using an asymptotic expansion 



cp' = + -ip 



1 



1 



,(2) 



r 0^ 

We get the approximation of the fixed points as 
ip* = mr + h~l) n {h -nit) + ..., 

P 



(8) 



o 
o 
o 
o 




FIG. 1: The left panel shows a schematic view of the electro- 
magnetic wave incident on the split ring Josephson junction 
resonator in the (y, z) plane. The right panel shows the equiv- 
alent circuit of the split-ring using the Resistively Shunted 
Junction model for the Josephson junction. 



where n is an integer. 

In the absence of damping a = and forcing h = 0, 
the system is Hamiltonian with 



1 



^,¥',r) = ^f T + /3(l 



x 1 2 1 2 
COS^) + -<p EE -p jT 



The stable fixed points correspond to the minima of the 
potential U(ip) and the even values of the integer n. Fig. 
[2] shows a plot of the potential U((f) for (3 = 1,9.76 and 
100. For j3 — 1 shown as a continuous curve (red on- 
line) there is only one fixed point ip — 0. For f3 — 9.76 
shown in dashed line (green online) there are three min- 
ima corresponding to stable fixed points, tp — 0, ±</>* 
where ip* « 2tt. For (3 — 100 there are many stable 
fixed points. 

For this one degree of freedom Hamiltonian, the orbits 
are the contour levels of the Hamiltonian. An important 
orbit is the separatrix connecting the two unstable fixed 
points up = ip* w 7r. It is given by LK = l K(ip*,Q). Fig. 
|3] shows the phase portrait for (3 = 9.76. For this value 
there are only five fixed points. Notice the closed orbits 
around the fixed points, the closed orbits surrounding the 
three stable fixed points. 

Another point is that the incident flux can be used to 
modify the energy levels of the system. Assuming the 
incident flux to be constant we can add a term to the 
potential and obtain the generalized potential 

U(ip) = 13(1 - cos <p) + X -ip 2 - hip, (10) 

where h is the normalized incident flux, assumed con- 
stant. This expression is plotted in Fig. 2] for f3 = 15 
and h = 0, 1.8-7T and 4.57T. The minima are symmetric 
for h = and they are shifted to the left and the corre- 
sponding value of the potential is decreased. By applying 
a sufficiently large continuous field one can then shift the 
system from one state to the other. Notice that for the 
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FIG. 2: Potential energy U(<p) = /3(1 - cos^) + \tp 2 with 
three different values of /?, /3 = 1, 9.76 and 100. 
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we will select the incident flux h to move the system from 
one equilibrium to another. 



III. SWITCHING BETWEEN EQUILIBRIUM 
STATES 



We consider now that the state of the system can be 
shifted from one fixed point to another via an incident 
flux. For a short lived perturbation, the system then 
relaxes freely to a minimum of energy. The influence of 
the damping is essential , it should be present to allow 
the relaxation but small to preserve the picture of the 
potential. To examine how an incident flux will shift 
the system from one equilibrium position to another it 
is useful to analyze the work equation. To obtain it, we 
multiply Q by ip^ T and integrate over time. We get the 
difference in energy 
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FIG. 3: Phase portrait (ip,ip T ) of the Hamiltonian system ([9]) 
(a — 0) for /3 = 9.76 and M — 0. The contour levels presented 
are 0.1, 1, 5, 17., 25.0360499136927 (separatrix) , 30. ,40. and 
50. 



large field h = 4.57T, there is a fourth minimum, around 
7tt. 

To summarize, the SRR-JJ has energy states U(ip) 
near 2nir. The number of these states is controlled by 
the parameter (3 i.e. the ratio of the flux to the flux 
quantum and the magnetic field h. In the next section, 
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FIG. 4: Potential energy U ((f) = /3(1 — cos <p) + ±ip 2 — hip for 
three different values of the static incident flux h — 0, 1.87T 
and 4.57T. 



E( T2 )-E(n) 



l 

2^ 



/3(l-cos(^) + ^ 2 



drhifir — a I drip 



(11) 



The first term on the right hand side is the forcing while 
the second one is the damping term. When a square 
pulse is applied to the system, such that 

/i(t) = a, for ri < r < T2, elsewhere 

the first integral is a\ip(r2) — ip(ri)]. If the system is 
started at (0, 0) in phase space so that the E(t±) = and 
</?(ti) = 0, we have 

E(t 2 ) = a<p(T 2 ) - a dripl, 

Jti 

so that f(T 2 ) determines how much energy is fed into 
the system. When the pulse is long t 2 » n <p(t) will 
relax and oscillate so that there are values of t 2 such that 
<p(t 2 ) 0. In that case no energy gets fed into the system. 
A sure way to avoid this is to take a narrow pulse. 

The natural frequency of the oscillator around the 
(0, 0) fixed point is 



w 



v^+l, 



(12) 



which for j3 = 9.76 gives ujq « 3.28 and a period 
Tb = 2tt/ujo Ri 1.91. To simplify matters we now con- 
sider a pulse of duration much smaller than To. This 
is experimentally feasible. For the present analysis the 
pulse duration is very small so that it can be approxi- 
mated by the Dirac delta function h(r) = aS(r), where a 
is a parameter. Let us analyze briefly the solution. The 
equation ((4]) becomes 



¥>tt + (a + S)ip T + (3sin(ip) + <p = aS(r). 



(13) 
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FIG. 5: Parameter plane (a, a) showing the different final 
states reached by the system, the left focus * (blue online), 
the center focus + (red online) and the right focus x (green 
online). The parameter /3 = 9.76. 



Integrating the equation on a small interval of size e 
around 0, we get 



[w]i e 



a\<P\- f + 



dr(P sin tp + ip) = a. 



(14) 



We now take the limit e — ► 0. We will assume continuity 
of the phase so that [p] e _ e -4 0. The third term being 
the integral of a continuous function tends to when 
the bounds tend to 0. Assuming <p T (0_) = we get 
yv(0+) = a so that such a short incident pulse will just 
give momentum to the system. With these initial condi- 
tions we solve equation Q numerically using a Runge- 
Kutta algorithm with step correction of order 4 and 5. 

We will now explore systematically the plane {a, a) 
characteristic of the incident pulse. The plot in the 
(a, a) parameter plane shown in Fig. [5] shows the fi- 
nal states, O the central focus (+), R the right focus (x) 
and L the left focus * reached by the system. Notice 
how these are organized in "tongues" following the se- 
quence OROLOROL... as one sweeps the plane 
counter clock- wise starting from the horizontal axis. This 
simple geometrical picture can be understood by examin- 
ing Fig. [31 In the case of small damping, the separatrices 
around the fixed points are not affected very much. Their 
perimeter is proportional to the probability of reaching 
one fixed point or another. The system is moving clock- 
wise along the orbits. Assume the system reaches the 



central point O for a given set of parameters. If the 
damping is increased, the orbit might not reach O but 
will settle in L. Similarly if more kinetic energy is given 
to the oscillator, it might reach R instead of O. 

Another important point is that the impulse given to 
the resonator must be very short so that it relaxes fol- 
lowing a free dynamics. The typical frequencies of these 
devices are about 500 GHz so the impulse must be around 
5 Thz. 

Now that we have seen how to prepare the system in a 
given state, it is important to check if this state is really 
reached. For that one can use a small signal to analyze 
the state by reflection or transmission. This is the object 
of the next section. 



IV. MICROWAVE SPECTROSCOPY OF THE 
SRR-JJ 

We assume that the ring is submitted to a fixed mag- 
netic field h s and that it is in a local minimum (p s of the 
potential. Then we send in a small electromagnetic pulse 
SH and examine the response Sip of the ring using the 
scattering theory. The linearized equations for Sh, Sip 
read 



Sh cc - Sh TT = -7 (-S^ TT + Sh tt ) 5(0, 
S$ TT + aS$ T + (3 cos(p s )5§ 4 8$ = Sh. 

We now assume periodic solutions 

Sh=h a e iuT , &<p = <p 9 e iuT , 

and obtain the reduced system 

haXC + u2h a = l^ 2 (-fa + h a ) S{Q, 

[-w 2 + iauj + 14/9 cos(ip s )] <p a = h a . 



(15) 



(16) 



(17) 



In the scattering we assume the electromagnetic wave to 
be incident from the left of the film located at f = 0. We 
then have 

h a = e-^+rRe* 1 * ,C < ; h a = 7e~ iu<: ,C> , (18) 

where H is the amplitude of the reflected wave and T 
the amplitude of the transmitted wave. We have the 
following interface conditions at £ = 

ha(0~) = h a (0+), [h a , ( ]°t = (-¥>„(0) 4 MO)) • 

(19) 

They imply two equations for 51 and T from which we 
obtain the transmission coefficient 



7 : 



2(— to 2 4 14/3 cos cp s 4 iua) 



D 



the reflection coefficient 



-aju 2 4 ijuj(/3 cosp s — ui 2 ) 
D : 



(20) 



(21) 
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FIG. 6: Square of the modulus of the reflection coefficient |i?| 2 
as a function of the frequency u for j3 = 9.76, h = 0, ip s — 
for a large coupling 7 = 1. The dashed line (red online) 
presents the reflection coefficient |i?| 2 for a polarized dielectric 
film for comparison. 

and where the denominator is 

D = 2(-uj 2 +l+f3cosLp s )+a'-fUj 2 +i{2auj-juj((3cos(p s -uj' 2 )}. 

(22) 

The square of the modulus of R is 

m \d\2 ■ ( 23 ) 

As seen in section 3, the split-ring oscillator has a finite 
number of equilibria tp s depending on the parameter /3. 
We will assume that a — throughout the section. This 
damping does not change the width of the resonances, it 
only affects the bounds of \R\. For a = 0, < \R\ < 1 
while for a > 0, a < \R\ < 1 - a. When a = \3l\ 2 is an 
even function of oj, therefore we will only consider uj > 0. 

As a first example we consider h = and /3 = 9.76 
corresponding to the potential plotted in Fig. [2] The 
reflection coefficient is plotted in Fig. [5] for 7 = 1 and 
(f s = 0. Notice how the resonance is asymetric around a 
minimum u> s . The system is passing for low frequencies 
only. This line shape is a typical Fano resonance [16[ due 
to the nonlocal coupling associated to the ip TT term in 
the right hand side of the wave equation (|3|). To see this, 
we have reported in Fig. [5] the reflection coeffcient for a 
simple polarized dielectric film [ItJ ■ This simpler system 
is passing for all frequencies except for a narrow band 
around the resonant frequency lo s . The expression for the 
resonant frequencies uj s can be obtained by considering 
the minima of |3i| 2 . These correspond to the second term 
in the numerator of ([23]) being zero. We get 

uis = a/8 cos ip s . (24) 

Let us now examine the influence of the different pa- 
rameters on the reflection coefficient. Fig. [7] shows the 
reflection coefficient |i?| 2 for the two steady states <p s and 
7 = 1 (left panel) and 7 = 0.2 (right panel). For the large 
coupling the resonance is sharp although for large to we 



FIG. 7: Square of the modulus of the reflection coefficient 
|7?| 2 as a function of the frequency ui for j3 = 9.76 and H = 
for a large coupling 7=1 (left panel) and a small coupling 
7 = 0.2 (right panel). 
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TABLE I: Spectroscopy data (uj s , ip a ) for a split-ring oscillator 
with three steady states. The parameters are /3 = 9.76, h = 
1.8tt, 7=1, a = 0. 



recover \R\ = 1. The value of the steady state ip s also in- 
fluences \R\. In Fig. [7J the continuous curve corresponds 
to ip s = and the dashed line corresponds to ip s — ±5.66, 
the two additional local minima. The two curves are very 
close and the reflection coefficient does not indicate if we 
are in the left or right minimum. Notice also that the 
system only passes low frequencies for 7 = 1. On the 
contrary for 7 = 0.2 it is globally passing for all frequen- 
cies except for a narrow band around u> s . These plots 
would enable to extract the coupling parameter from a 
SRR-JJ. 

As a second example we consider /3 = 9.76 and h = 
1.87T corresponding to the potential plotted in Fig. [4] 
There are three minima given in the table I. The reflec- 
tion coefficient |i?| 2 is plotted in Fig. [5J Now it is possi- 
ble to distinguish the "left" minimum 1 from the "right" 
minimum 3 because the non zero magnetic field h has 
lifted the degeneracy. 

Finally we consider h = and (3 — 30 for which the 
potential U (4>) has five minima indicated in the table II. 
The square of the modulus of the reflection coefficient 
([2"B"f is plotted in Fig. [§] for the five different equilibria. 
In the example shown, the spectroscopy data (oj s ,ip s ) is 
given in table II. Again there is no possibility to dis- 
tinguish whether the system is in the left or the right 
minimum. 
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FIG. 8: Square of the modulus of the reflection coefficient 
\R\ 2 as a function of the frequency uj for the three dif- 
ferent equilibria indicated in Table I. The parameters are 
/3 = 9.76, h = 1.87T, 7 = I and a = 0. 



index 


1 


2 


3 


4 


5 







6.08 


12.15 


18.2 


24.18 




5.477 


5.420 


5.238 


4.888 


4.169 



TABLE II: Spectroscopy data (cj s , tp s ) for a split-ring oscilla- 
tor with five steady states. The parameters are j3 — 30, 7 = 
1, a = 0.01. 



system from the ground state to one of these excited 
states. By sending a microwave field on the resonator 
we can perform a spectroscopy of it and characterize in 
which state it is. Using a scattering theory formalism we 
computed the reflection and transmission coefficients for 
the wave. These have a specific minimum that we can 
compute analytically. Another interesting feature is the 
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FIG. 9: Square of the modulus of the reflection coefficient 
\R\ 2 as a function of the frequency ui for the five different 
equilibria indicated in Table II. The index of the minima is 
indicated in the plot. The parameters are /3 = 30, 7 = 1 and 
a = 0.01. 



V. CONCLUSION 

We considered the interaction of an electromagnetic 
wave with a split-ring Josephson resonator. For appropri- 
ately chosen parameters, the resonator has excited states 
whose number can be controlled by carefully tuning the 
inductance and capacity of the ring. The existence of 
these excited states makes this system similar to an ar- 
tificial atom with discrete energy levels. Other artificial 
atoms containing nonlinear elements like a diode Q, a 
Kerr material U or a laser amplifier [6] would not give 
these discrete levels. In addition, since the oscillator is 
operating in the superconducting regime the losses are 
very small as opposed to the current metamaterials. 

First we assumed that there are just two excited states 
and showed how an incident magnetic flux can shift the 



form of the resonance which is Fano like and very differ- 
ent from the one obtained for a ferroelectric thin film [ljj ■ 
Note that the right and left minima are undistinguishable 
with this measurement and that a static magnetic field 
will lift this degeneracy. We hope this analysis will be 
useful to experimentalists studying split ring Josephson 
resonators. 
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